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ABSTRACT 

We examine the effects of significant electron anti-neutrino fluxes on hydrogen 
burning. Specifically, we find that the bottleneck weak nuclear reactions in the 
traditional pp-chain and the hot CNO cycle can be accelerated by anti-neutrino 
capture, increasing the energy generation rate. We also discuss how anti-neutrino 
capture reactions can alter the conditions for break out into the rp-process. We 
speculate on the impact of these considerations for the evolution and dynamics 
of collapsing very- and super- massive compact objects. 

Subject headings: neutrinos — nuclear reactions, nucleosynthesis, abundances - 
stars: evolution 



1. Introduction 

Hydrogen burning involves the conversion of four protons into an alpha particle, two 
positrons, neutrinos and photons. The principal bottleneck involved in this process is the 
weak interaction conversion of protons into neutrons. For decades the primary mechanisms 
of hydrogen burning have been an astronomical staple. Bethe & Critchfield (1938) first 
elucidated the proton-proton chain (pp-chain) where the weak conversion is accomplished 
by two protons interacting to become a deuteron, p(jp, v e e + )d. von Weizsacker (1937, 1938) 
and Bethe (1939) independently described the CNO cycle, where carbon is used as a catalyst 
in hydrogen burning, and the weak conversion of protons to neutrons occurs through the 
positron decay of isotopes of oxygen with half lives of about 100 seconds. 

A large flux of electron anti-neutrinos (u e ) could alter the hydrogen burning paradigm. 
Anti-neutrino capture could perform the necessary conversion of protons to neutrons. The 
^-capture cross sections of relevance are very small, but depend strongly on neutrino energy. 
The smallness of these cross sections allows energetic neutrinos to escape from deep within a 
compact object, where the temperature and other energy scales are high, and freely stream to 
where hydrogen burning is occurring. Nevertheless, if i/ e -capture is to have a significant effect 
on hot hydrogen burning, a truly prodigious flux (0 Pe > 10 40 cm -2 s -1 ) and large neutrino 
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energy ((E Pe ) > a few MeV) would be necessary. It should be kept in mind, however, that to 
affect hydrogen burning, the z/ e -capture rates need only be comparable to the corresponding 
positron decay rates. 

The difficulty would be to find an environment capable of producing these fluxes of 
neutrinos, yet quiescent enough that simple hydrogen burning could be relevant, and the 
products of such burning could be ejected into space. High entropy electron-positron plasmas 
are efficient engines for the production of neutrinos and anti-neutrinos of all flavors. Possible 
environments that may merit future investigations into the effects of anti-neutrino capture on 
hydrogen burning include high mass accretion disks and collapsing very- and super- massive 
objects. 

In this paper, we investigate the effects of a prodigious neutrino flux on hot hydrogen 
burning. In section 2 we point out the effects of anti-neutrino capture on the rate limiting 
steps in both the pp-chain and the /5-limited CNO cycle, and its implications for the energy 
generation rates. In section 3, we examine the consequences for the rp-process and energy 
generation mechanisms. In section 4, we consider the case of a supermassive star collapsing 
on the general relativistic Feynman-Chandrasekhar instability, and the effects of its internal 
neutrino production on hydrogen burning in its envelope. We give conclusions in section 5. 

2. Neutrino-Induced Hydrogen Burning Mechanisms 

The rate limiting step in the pp-chain is the weak interaction conversion of two protons 
into a deuteron, a positron, and an electron neutrino. A significant flux of electron anti- 
neutrinos allows an alternate mechanism to be favored, where anti-neutrino capture on a 
proton creates a neutron and a positron {p e +P — > n+e + has been considered in supermassive 
objects by Woosley (1977) and Fuller & Shi (1997)). This step would be followed by a fast 
radiative proton capture to form a deuteron. Comparing the two reaction rates: p(p, v e e + )d 
vs. p(u e ,e + )n(p, , ~f)d, we find that for the prodigious anti-neutrino fluxes discussed in the 
introduction (0 Pe > 10 40 cm -2 s _1 , (E Pe ) > a few MeV) the anti-neutrino capture path is 
significantly faster in relevant astrophysical environments. This provides not only a new 
reaction path for hydrogen burning, but increases the energy generation rate by several 
orders of magnitude. 

The /5-limited CNO cycle, or hot CNO cycle, proceeds at a rate dictated by the positron 
decay of 14 and 15 0, with half lives of 71 s and 122 s, respectively. (See for example Hoyle 
& Fowler 1965 and Audouze, Truran, & Zimmerman 1973.) These decays likewise could be 
augmented by electron anti-neutrino capture, 14 0(z/ e , e + ) 14 N and 15 0(z/ e , e + ) 15 N. Figure 1 
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shows the acceleration of the relevant weak rates as a function of total electron anti-neutrino 
flux for an assumed Fermi-Dirac z/ e -energy spectrum with average z/ e -energy (E Pe ) = 10 MeV 
and zero chemical potential. The flux at which anti-neutrino capture becomes important 
scales appoximately as (E Pe )~ 2 . For a large enough flux (0 Pe > 10 39 cm -2 s^ 1 in the case of 
Figure 1), the reaction rates are proportional to the incident flux of electron anti-neutrinos. 
Our weak rate calculations are described in Appendix A. 

Additionally the CNO cycle is accelerated by the presence of free neutrons. The strong 
interaction reactions 15 0(n, p) 15 N and 14 0(n,p) 14 N have a significantly larger cross section 
than the electromagnetic reaction n(p,'~f)d. As a result, neutrons are diverted from the 
modified pp-chain into the CNO cycle. Figure 2 shows how the neutrons created by pipe-, e + )n 
are distributed between the competing reactions 15 0(n,p) 15 N, 14 0(n,p) 14 N and n(p^)d. 
Notice that for an assumed Fermi-Dirac z/ e -energy spectrum with (E Pe ) = 10 MeV and 
zero chemical potential, the ratio of neutron captures on 15 to 14 to p is approximately 
4.5 : 2 : 1 for a large range of z/ e -fluxes. 

Figure 3 illustrates the most significant reaction flow paths involved in hydrogen burning 
when a significant z/ e -flux is present. The pp-chain is modified as anti- neutrino capture allows 
the circumvention of the slowp(p, u e e + )d reaction. Also included are the triple-alpha process, 
which would provide a path between the pp-chain and CNO cycle, and the break out into 
the rp-process via 15 0(o;, 7) 19 Ne(p, 7) 20 Na. (Wallace & Woosley 1981) 



3. Side Effects 

A large flux of electron anti-neutrinos certainly accelerates the weak rates that provide 
the bottleneck in hot hydrogen burning. However, since this flux also increases the rates of 
other positron decays, a number of side effects are possible. 

A principal mechanism for break out into the rp-process involves the reaction path 
15 0(a, 7) 19 Ne(p, 7) 20 Na. The criteria for break out into the rp-process can be found in 
the competition between proton capture on 19 Ne, and the decay of 19 Ne through positron 
emission and now, anti-neutrino capture. Thus for densities and temperatures that satisfy 
the inequality 

pXA P7 ( 19 Ne) > A w ( 19 Ne), (1) 

break out into the rp-process will occur (Wallace & Woosley 1981). Here the density p is 
in gcm~ 3 , X is the hydrogen mass fraction, Au,( 19 Ne) is the total weak decay rate of 19 Ne 
(positron emission and z/ e -capture), and A P7 = Na(cv) p1 , where Na is Avagadro's number 
and the thermally averaged product of cross section and speed is taken from Caughlan & 
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Fowler (1988). 

Including a large flux of electron anti-neutrinos would result in higher weak decay rates 
(X w ). This increase in the right hand side of equation (1) would require an increase in 
temperature (increasing A P7 ( 19 Ne)) for a given density at which break out into the rp-process 
would occur. Figure 4 shows the effects of an electron anti-neutrino flux on the conditions 
necessary for break out into the rp-process. 

The pp-chain is the dominant process of energy generation in the sun, while the CNO 
cycle is dominant in stars that are more massive. However, with a large flux of electron anti- 
neutrinos, these processes become independent of temperature so long as the temperature 
is high enough to guarantee that proton capture remains comparatively fast. A high flux 
of electron anti-neutrinos allows the pp-chain to compete favorably with the CNO cycle. 
For example, Pruet et al. (2005) have studied nucleosynthesis in supernova winds where 
hydrogen "burning" is completely dominated by u e and v e capture on free nucleons (see also 
Qian et al. 1993). 

The scarcity of 15 in comparison to free protons means that for large anti-neutrino 
fluxes and average energies, the pp-chain is the dominant mechanism in hydrogen burning at 
temperatures that the CNO cycle would traditionally dominate. Figure 5 shows a comparison 
between the energy generation rates of the pp-chain and the CNO cycle for X/Z' = 10 and 
100, where X is the hydrogen mass fraction and Z' is the mass fraction in carbon, nitrogen, 
and oxygen isotopes. For large anti-neutrino fluxes and average energies, the pp-chain is the 
dominant energy generation mechanism, while for low fluxes and average energies the CNO 
cycle takes over. 



4. Example: Supermassive Stars 

Now we consider the case of a supermassive star, a star so massive that it collapses on the 
general relativistic Feynman-Chandrasekhar instability (M >5x 10 4 M Q ) (Fuller, Woosley, 
& Weaver 1986, hereafter FWW). If such objects did exist, for example in the early universe, 
their homologous cores would emit copious fluxes of neutrinos and anti-neutrinos of all flavors 
during their collapse (Shi & Fuller 1998). Shi and Fuller examined the collapsing core of 
a supermassive star, calculating the luminosity and energy spectrum of neutrinos emitted. 
The total neutrino luminosity was found to be 

L v « 2.8 x 10 57 {Mf c Y l - b erg s"\ (2) 

where Mf c is the mass of the homologous core in units of 10 5 M Q . Additionally they found 
the energy spectrum of neutrinos of all flavors to fit remarkably well to a Fermi-Dirac spec- 
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trum with a higher temperature than the central plasma temperature (T„ ~ 1.6T) and a 
degeneracy parameter (chemical potential divided by temperature) r\ v fa 2. 

We can check the effect of this flux of neutrinos an antineutrinos on the nuclear physics 
in the gas in the envelope of the star. As a point of reference, we choose a radius of 100 
Schwarzschild radii (r = 3 x 10 12 M^ C cm) where the gravitational binding energy of a 
nucleon is approximately equal to the nuclear energy liberated in these reactions, so there 
is a chance that any new nuclear physics that occurs as a result of z/ e -capture could be 
relevant. By "relevant" we mean that it is conceivable that material from this location could 
avoid being swallowed by the black hole forming in the core. Only a detailed simulation 
with general relativistic hydrodynamics could reveal whether or not material affected by 
z/ e -capture is ever ejected into space. We are now free to repeat the analyses done above, 
but with one free parameter, the mass of the collapsing homologous core. 

Figure 6 shows the acceleration of the relevant weak rates as a function of homologous 
core mass. The anti-neutrino capture rate is proportional to (Mf 10 ) -4 . We see that the effects 
of anti-neutrino capture on the decay rates of 14 and 15 become insignificant for large 
homologous core masses {Mf c > 0.4). If supermassive stars ever formed, it is conceivable 
that they were in the first generation of stars with primordial initial abundances. In this case 
the CNO cycle would be negligible, at least initially (FWW). However, even in this case, the 
energy generation rate of the pp-chain would be boosted by several orders of magnitude. It 
would be interesting to see if this added energy source would have a discernible effect on the 
eventual fate of a collapsing supermassive star. 

Including a hydrogen burning phase in the final stages of the collapse of a supermassive 
star may affect the eventual fate of its baryons. Hydrodynamic, post-Newtonian calculations 
done in FWW show that initially metal-free supermassive stars will collapse to black holes. 
Shibata & Shapiro (2002) use a fully relativistic simulation in axial symmetry to deduce 
that the supermassive star collapses to a black hole surrounded by some remaining gas in 
an ambient disk. 

The principal nucleosynthetic issue is whether any material that had experienced v e 
capture-affected hydrogen burning escapes being incorporated into a black hole. Of course, 
there is the prior issue of whether material at the relatively low temperatures and densities 
which characterize hydrogen burning ever experiences high z/ e -fluxes. Both issues are related: 
to see nucleosynthesis products of v e capture-affected hydrogen burning, the material must 
be ejected before the point at which nuclear burning proceeds past simple hydrogen burning 
and approaches or attains nuclear statistical equilibrium (NSE). We are skeptical that these 
conditions can be met. Fully relativistic numerical simulations could settle these issues. 
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Obviously, the NSE nucleosynthetic yield is uninteresting in the context of this paper. 
However, a mass shedding scenario could be conducive to conditions that favor hydrogen 
burning and the rp-process. Speeding up weak decays could affect the relative abundances of 
the rp-process elements. A simulation that follows these species and their chemical reactions 
would be necessary to address this issue. 

5. Conclusions 

In this paper we have examined the effects of a prodigious flux of electron anti-neutrinos 
on hydrogen burning. We have found that the traditional positron decay bottlenecks in 
hydrogen burning can be removed and replaced by much faster z/ e -capture reactions under 
some conditions. This would result in an increase of several orders of magnitude in the 
energy generation rate over what would be expected without such a flux. 

Additionally, the v e -Qux would alter the conditions necessary for break-out into the re- 
process, increasing the temperature necessary to do so at a given density. If conditions allow 
the break-out into the rp-process, we could expect an acceleration of the flow toward the 
iron-peak facilitated by and accelerated by z/ e -capture. 

When applied to the neutrino flux emitted in the final stages of the collapse of a su- 
permassive star, interesting changes from current simulations may occur on the lower end of 
the supermassive star mass spectrum. Whether or not these effects are relevant, remains an 
open question that can only be answered by simulations that are able to include hydrogen 
burning during the final collapse of the star. 

Important issues that remain open include finding an astrophysical environment where 
the effects discussed here could take place. Accretion disks surrounding black holes may 
provide a combination of high accretion rates and hot, high entropy disks which could pro- 
duce the necessary fluxes of electron anti-neutrinos. (See for example Surman, McLaughlin, 
& Hix (2005) for a discussion of neutrino emission in lower mass accretion disks.) Super- 
massive stars may exhibit this effect, though there is uncertainty related to whether or not 
these objects ever existed. Computer simulations would be useful to determine any changes 
in the expected nucleosynthetic yield, and the effects of their possible distribution into the 
surrounding IGM. 

We would like to thank S. E. Woosley, Y.-Z. Qian and A. Heger for useful discussions. 
This work was supported in part by NSF grant PHY-04-00359 and the TSI collaboration's 
DOE SciDAC grant at UCSD. C.T.K. would like to acknowledge a fellowship from the ARCS 
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A. Calculation of Weak Rates 

In this work we calculate the z/ e -capture rates in the manner described in Fuller, Fowler, 
& Newman (1980, 1982) (FFN I, II) and Fuller & Meyer (1995). We employ measured 
discrete states only. 

Our 14 0(z/ e , e + ) 14 N rate calculation includes only the 14 ground state (spin and parity 
J w = + ) and the measured weak branches to the 14 N ground state (J 71- = + , log 10 ft = 7.3), 
first excited state (J 71 " = 1 + , log 10 ft = 3.5), and second excited state (J 71 = 1 + , log 10 ft = 
3.1). Contributions to the stellar rate from thermal excitation of parent states are small here 
as a result of the high first excited state excitation energy (5.17 MeV) and the temperatures 
of interest. Likewise, branches to higher excited states in 14 N are not significant. A possible 
exception is the first isobaric analog state in 14 N (J n = + ) at excitation energy of 8.62 MeV. 
This branch will have a large matrix element but will be Q-value-hindered relative to the 
+ — > + ground state to first excited state, pure Fermi branch. 

Our 15 0(z/ e , e + ) 15 N rate calculation includes only the ground state (J n = l/2~) branch. 
This channel has a large matrix element, corresponding to log 10 ft = 3.6. Branches to 15 N 
excited states will not be significant. 15 N states below 9.15 MeV excitation energy have 
positive parity and the branches from the 15 O ground state to them will be forbidden. We 
note, however, that 15 and 15 N are isospin mirrors. This can be a significant fact for stellar 
weak interaction rates, as it implies large Fermi and Gamow- Teller matrix elements coupling 
each parent state with its daughter isobaric analog state (FFN I, II). Thermal excitation of 
15 excited states would open weak branches to corresponding isobaric analog states in 15 N. 
This is not likely at the temperatures of interest because the first excited state of 15 is at 
about 5.2 MeV excitation. 

Our calculation of the 19 Ne(z/ e , e + ) 20 Na rate includes only the ground state of 19 Ne 
(J 7 " = 1/2+) and branches to the ground (J 7 " = 1/2+) and third excited state (J w = 3/2 + ) of 
19 F. The first of these branches, with log 10 ft = 3.2, dominates the rate. We note, however, 
that 19 Ne and 19 F are isospin mirrors. Since temperatures are high near CNO cycle breakout, 
thermal excitation of the first (J n = 5/2 + ) and second (J n = l/2~) excited states can be 
expected to carry a fraction of the total weak rate. However, on the assumption that the 
matrix elements for these branches are identical to that for the ground-to-ground transitions, 
inclusion of these branches makes little difference (< 1%) for the rates and our conclusions. 
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Fig. 1. — Key weak decay rates as a function of electron anti-neutrino flux, assuming a 
Fermi-Dirac z/ e -energy spectrum with zero chemical potential and average z/ e -energy (E Pe ) = 
10 MeV. The solid lines, labeled \ u and A i5 are the total weak decay rates of 14 and 15 0, 
respectively, through positron decay and anti-neutrino capture. The dashed line is the rate 
of conversion of a proton to a neutron via anti-neutrino capture alone. 
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Fig. 2. — Fraction of neutron captures on 5 (/15), O (/i 4 ) and p (f np ), assuming a Fermi- 
Dirac i> e -energy spectrum with zero chemical potential, average z/ e -energy (Ep e ) = 10 MeV 
and X/Z' = 100, where X is the hydrogen mass fraction and Z' is the mass fraction in 
carbon, nitrogen and oxygen isotopes. 
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Fig. 3. — Qualitative picture of hydrogen burning under the influence of a prodigious elec- 
tron anti-neutrino flux. Shown are a modified pp-chain, /3-limited CNO cycle, and the less 
significant triple-a reaction and break-out to the rp-process. 
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Fig. 4. — Conditions for break-out into the rp-process, assuming a Fermi-Dirac z/ e -energy 
spectrum, zero chemical potential, and (E Pe ) = 10 MeV. The plotted contours are for, 
in ascending order, <p Ve = 10 38,40,42 cm~ 2 s _1 . Zero neutrino flux is indistinguishable from 
(f)- e = 10 38 cm" 2 s -1 . 
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Fig. 5. — Comparison of energy generation rates for the pp-chain and CNO cycle at a 
temperature of 5 x 10 8 K. The contours are where the energy generation rates are equal 
for the two hydrogen burning mechanisms for values of X/Z' = 10, 100 (from right to left). 
Larger i/ e -fluxes and average energies favors the pp-chain, while smaller fluxes and average 
energies favor the CNO cycle. 
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Fig. 6. — Rates of relevant weak decays at a distance of 100 Schwartzschild radii from the 
center of a supermassive star with homologous core mass M^ c x 1O 5 M . As in Figure 1, the 
solid lines are the total rates of decay of 14,15 0, and the dashed line is the proton to neutron 
conversion rate through anti-neutrino capture alone. 



